Modified Newtonian Dynamics In Dimensionless Form 



in 
o 
o 

(N 

Oh- 
< 



o 
o 



Oh- 



X 



W. F. Kao* 

Institute of Physics, Chiao Tung University, Hsinchu, Taiwan 

Modified Newtonian dynamics proposed that gravitational field needs modifications when the 
field strength g is weaker than a critical value go- This has been shown to be a good candidate 
as an alternative to cosmic dark matter. There is another way to look at this theory as a length 
scale dependent theory. One will show that modification of the Newtonian field strength depends 
on the mass distribution and the coordinate scale of the system. It is useful to separate the effective 
gravitation field g(r) into a small scale (or short-distance )g s field and a large scale (or a long- 
distance) gi field that should be helpful for a better understanding of the underlying physics. The 
effective potential is also derived. 



I. INTRODUCTION 



Modified Newtonian dynamics (MOND) was proposed by Milgrom (1983) that gravitational field needs modifications 
when the field strength g is weaker than a critical value go. This has been shown to be a good candidate as an 
alternative to cosmic dark matter (Sanders 2001). The phenomenological foundations for MOND are based on two 
observational facts: (1) flat asymptotic rotation curve, (2) the successful Tully-Fisher (TF) law, M ~ V a (Tully & 
' Fisher 1977) for the relation between rotation velocity and luminosity in spiral galaxies. Here a is close to 4. 
O'N In this paper, one proposes a different view of MOND by looking at the physics related by the mass distribution 

and the coordinate scale of the system. In addition, one finds it useful to separate the effective gravitation field g(r) 
into a small scale (or short-distance )g s field and a large scale (or a long-distance) gi field that should be helpful for a 
better understanding of the underlying physics. The effective potential is obtained. Possible relation with the induced 
| gravity model is also discussed. 

If dark matter does not exit, rotation curve observations indicate that g behaves as 1/r, a typical 2-dimensional 
attraction field, at large distance greater than the galactic scale around 10 5 lightyears. At short distance, g goes 
like 1/r 2 , a typical behavior of a 3-dimensional attraction field. Therefore, it is also very interesting to study the 
changing pattern of the effective dimension d(r) which will be defined as a function of the physical scale r. One will 
plot this d(r) functions hoping for a better understanding of the changing pattern of the effective dimension. Possible 
speculation with the Kaluza-Klein theory will be discussed in this paper too. 

We will start by showing how to simplify the expression of MOND field strength by writing it with suitable scale 
parameters. Our aim is to write all physical quantities in dimensionless form with suitable physical units. 

First of all, it was pointed out (Milgrom 1983) that there exists a critical acceleration parameter go = 1.2 x 
10 _8 cms -2 characterizing the turning point of the effective power law associated with the gravitational field in 
■ MOND. Gravitational field of the following form was suggested 

a | 

9-^{ — )=9n (1) 
5o 

with a function /i considered as a modified inertial. Here g^ is the Newtonian gravitational field produced by certain 
mass distribution. Milgrom argues that 

<*> = 7^ < 2 > 

provides a best fit with many existing observation data including the rotational curve of many spiral galaxies. Ob- 
servations are selected mostly from the most reliable 21 cm hydrogen line measurements. The 21 cm hydrogen line 
is the photon emission due to the superfine energy difference resulting from the magnetic dipole-dipole interactions 
between the cored proton and orbiting electron. The existing neutron hydrogen all around the invisible region of 
spiral galaxies provides better estimates on the rotation curves beyond the visible region of spiral galaxies. 
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Milgrom also points out that the critical parameter go is close to the value of cHo ~ 5 x 10 _8 ctos -2 . It is also 
the gravitational field produced by an to — 200MeV massive particle at its Compton surface. Indeed, one finds that 
g200MeV = Gm/X = Gm 3 c 2 /h 2 <~ 1.2 x I0~ 8 ctos~ 2 . Milgrom hence speculates that MOND could have to do with 
the quantum theory of the universe. 

Milgrom's original idea is that the gravitational field needs modification when g is weaker than critical field strength 
.go- It is interesting to find that r ~ 49700Zy if we define the scale length r via the following equation 

GM 

9o = (3) 
r o 

with M = 2 x 10 n M s and M s the solar mass. Note that one takes Mq close to the order of magnitude of a typical 
galaxy similar to our Milky Way. Taking a different value of M will only change the value of r . Selecting different 
value of M and hence ro will not affect our arguments in this paper. One chooses M and hence the corresponding 
ro simply to set them as units of scale. Physics will not be affected by the choice of unit. 

One observes, however, that r is pretty close to the radial size of the visible boundary of our Milky Way if M is 
chosen close to the total visible mass of Milky Way. This is certainly true as modification is needed only at large scale 
according to MOND. Starting next section, one will try to write the modified field strength g in a dimcnsionless form 
with the help of the parameters Mq, tq and go. The idea is to write M, r and g in units of Mo, ro and go respectively. 

II. MODIFIED NEWTONIAN DYNAMICS 

Milgrom argues that one can take the effective inertial [i as /i(x) = xj\J\ + x 2 . Therefore, one obtains 

by dividing both sides of the equation [I] by the critical parameter g . One can therefore write above equation as 

9N (5) 



y/T+9- 

with g and gN written in unit of go. This makes g and gN dimcnsionless from now on. 

Throughout this paper, we will focuss on the study of the system with a Newtonian attraction of the form gN = 
Gm/r 2 . This is the field strength at a radial distance r from a spherical distributed matter with total mass to. System 
with different mass distribution can be obtained straightforwardly. Note that Gm/(r 2 go) — mr^/{Mor 2 ). Therefore, 
one can write gN — ra/r 2 if gN, r and to are written in unit of go, ro and Mq respectively. Therefore, the modified 
field strength g can be written as 



m 

„2 



(6) 



One can further suppress the parameter to by defining r c = ro\fm and write r as a dimensionless coordinate variable 
in unit of r c instead of r . As a result, one can write above equation in a very compact form 

» 3 (7) 



in terms of above dimensionless parameters. To emphasize again, one writes 

TO 

tt > m, 

Mo 

9 



► 9, 

90 



r 



r c Jmro 



(8) 



in order to simplify the expression of the modified equation and make g, gN, and r dimensionless for convenience. It 
is straightforward to restore all dimension parameters when one needs to evaluate any corresponding physical values. 
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First of all, Eq. (7) can be solved to write g(r) as a function of r: 



Vl + Vl + 4r- 4 
y/2r 2 



= -exp[ 



sinh -l^„2^-l 



L (2r 2 



(9) 



One can plot the functions g(r) and compare it with the functions 1/r and 1/r 2 for a better picture of its small scale 
and large scale behaviors. Indeed, Figure 1 shows the behavior of g(r), 1/r, and 1/r 2 together for comparison. Because 
the difference at short distance scale is too small for comparison, one also plot it in logarithmic scale. Figure 2 was 
plotted as logarithm of these functions to enlarge the small difference among these functions in the short distance 
region. 
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FIG. 1. g(r),l/r and 1/r 2 plotted as thick solid line, dashed line, and thin solid line respectively. Note that g and r are 
dimensionless quantities defined in Eq. [8] . 
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FIG. 2. Logarithm of g(r),l/r and 1/r 2 plotted as thick solid line, dashed line, and thin solid line respectively. Note that g 
and r are dimensionless quantities defined in Eq. [8]. 

It is apparent that g(r) goes like 1/r 2 at short distance scale where r < 1. On the other hand, g(r) goes like 
1/r at large distance scale where r»l. One can also easily tell from the figures that they agree with our analytic 
approximation at both the regions r <g; 1 and r» 1 one obtained earlier. 

More importantly, writing g(r) as a function of r — * r/(y / roro) clearly shows that the Newtonian gravitational 
field receives significant modification only when r > 1, or equivalently, r > r c = ^JmrQ. This is in fact equivalent to 
the assumption made by Milgrom that significant modification is required when g < go- The picture looking at the 
physics related to the physical scale is however easier for a better understanding of possible underlying physics. 

Indeed, one easily finds that, if Eq. (9) is the correct gravitational field for all scales, the effective dimension of the 
system seems to decrease from 3 to 2 as the scale r increases. Geometric dimension seem to decrease if one increase the 
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physical scale of interest. The major difference is that the critical scale depends not only on r$ = <J GMq /go but also 
proportional to y/m, the square-root of the total mass of the system. One recalls that Kaluza-Klein theory proposes 
that observable physical dimension of a system increases as the physical length scale decreases, or equivalently, the 
energy scale increases. One will be able to observe higher dimensional effect at shorter distance scale according to 
Kaluza-Klein's original idea. 

The situation here is quite similar. When the physical length scale increases, the effective physical dimension 
appears to decreases from 3 to 2. Moreover, the critical length scale r c = \/mro depends on the total mass of the 
system. If m increases, equivalent to the increasing of total energy of the system, MOND modification is only apparent 
at larger scale. Imagine one approaches a system from a large distance area. One will be able to see the effect of one 
more dimension as one approaches closer to the center of the system. Increasing m will make one easier to see the 
differences at the same physical length scale. This agrees with the central idea of Kaluza-Klein approach. One can 
possibly say that energy or mass opens up more dimensions if Eq. (9), or any similar expression of the MOND field 
strength, is the correct equation for the gravitational field. 

Note that one can also tell from Eq. 9 that one can expand the effective theory in r if one needs analytic form of 
the theory either at short or large scale. The expansion in r is in fact the expansion in r/r c . When the scale length 
r/r c is small, it returns to the conventional Newtonian model. It will receive significant modification when r is large 
compared to r c . Therefore, attention should be addressed to the physical origin of the parameter r c . 

In order to take a more close look at the changing pattern of the above g, one can split it into two different parts: 



9 = 9i+9s = 



2y/2r 2 



Vl + 4rVl + Vl + 4r 4 \/2rVl + 4r 4 ' 



(10) 



Note that the first term is gi that goes like 1/r when r 1, while the second term is g s that goes like 1/r 2 when 
r <C 1. Therefore it is easy to see that gi and g s represents the long distance and short distance field strength of the g 
respectively. One is hoping that successful separation of g May shed more light in the search of the underlying theory. 
In fact, one can show that 
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(ii) 

(12) 



and hence, 



:(1 + ^ 4 ) 



(13) 



as r — ► 0. This indicates that the first order correction to the short-range gravitational field g is 0(r 2 ) in comparison 
to the 1/r 2 Newtonian interaction. Hence the first order correction is fourth order in r. This correction is very tiny 
if one takes the mass m to be the solar mass M s . In this case the critical radius r c O.lllZy. Hence, the first order 
correction 8g at a distance from our sun near r = 10 AU, would be around a small difference near 4.16 x 10~ 12 x g w . 
Here gio denotes the Newtonian gravitational field at a distance r = 10 AU from our sun. This would probably large 
enough to be measured in the near future. 

In Figure 3 one plots g(r),gi(r), and g-z(r) for comparison. We also plot the logarithm of these functions in Fig. 4 
in order to signify the small difference in short distance. 
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FIG. 3. g(r),gi(r) and g s (r) plotted as thick solid line, dashed line, and thin solid line respectively. Note that g and r are 
dimensionless quantities defined in Eq. [8]. 
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FIG. 4. Logarithm of g{r),gi(r) and g a (r) plotted as thicker solid line, dashed line, and solid thin line respectively. Note 
that g and r are dimensionless quantities defined in Eq. [8]. 



In addition, one can also compute the changing rate of the gravitational field for the following equation of running 
gravitational field 



2 (.9 2 + l) 3/4 



If we define k = 1/r as a parameter simulating the effect of the conjugated momentum, one can also show that the 
changing rate dkg{k) obeys the following equation 



,( g 2 + l) 5/ 4 
5 2 + 2 



d k 9 = 2 ^ J 7 o (15) 



These equations play the role similar to the renormalization group equation for a quantum field theory. In addition, 
Eq. [14] also states that g(r) decreases monotonically as r increases. Similarly, Eq. [15] indicates that g(k) increases 
monotonically as k increases. 



III. EFFECTIVE POTENTIAL AND EFFECTIVE DIMENSION 



One can integrate g for the effective potential. After some algebra, One can show that the effective potential 
Vi = — / g; • dr — J r gidr and similarly for V s = f r g s dr can be evaluated directly to give 



^ = m^ i+ ^ +4r4 -f: ^ k -y (16) 

V2 ^ 2n ■ n\2 n (l + Vl + 4r 4 )" 



and 



One can verify directly that V/ = gi and — g s and prove that above equations are indeed correct up to an irrelevant 
integration constant. In fact, it is difficult to specify this constant of integration in the conventional approach which 
take V(r — > oo) — > 0. This is because the effective potential in fact diverges at spatial infinity due to the logarithm 
behavior of the dominating 2D-like potential. 
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FIG. 5. V(r),Vi(r) and V s (r) plotted as thick solid line, dashed line, and thin solid line respectively. 



In order to take a look at the form of these effective potential , one plots V, Vj, and V s together in Fig. 5. It is easy 
to find that Vi has a flat plateau area when r — ► indicating that <?; contributes little to g at short distance. Similarly, 
one finds that the slope of V s becomes constant as r ^> 1 indicating that g s is significant only at short distance scale. 
Note that, the series sum in V s converges quickly. One takes n = 100 when we plot Fig. 5. 

One knows that the four dimensional space time is modified when the physical scale becomes small in the Kaluza- 
Klein theory. Recent evidences in MOND research seems to imply that the effective dimension of our universe decreases 
as the length scale increases. It is possible that Kaluza-Klein approach could play an important, but quite different, 
role for large scale system. The observed rotation curve of spiral galaxies indicates that the gravitational attraction 
goes as 1/r in all directions without the existence of any cosmic dark matter. This is, however, quite different from 
the conventional geometry of any two dimensional space. But the success of MOND seems to indicate a new direction 
of dimensional effect. 

It is then very interesting to study how the effective dimension decreases as the observation scale increases. There- 
fore, one manages to define an effective dimension in this paper hoping to find our way to the discovery of the 
underlying theory. 

Indeed, one can define the effective dimension d(r) by the definition 



1 



n d-l 



(18) 



such that d — > 3 as r — > and d — > 2 as r — > oo. This is a naive way to define an effective dimension such that d will 
signify the geometric dimension as 3-dimension and 2-dimension respectively in a close way. It is then straightforward 
to show that 



-d—3 



y/l + Vl + 4r 4 

71 



and hence, 



d= 2- 



sinh 



2r 2 



21nr 



(19) 



(20) 



Note, however, that this definition of d(r) fails to tell us anything when r — ► 1. The reason is that d is ill-defined at 
r = 1 which is an obvious result of Eq. 18. One resolution to this problem is to introduce a A-shaped smooth function 

A(r) that peaks at r = 1 reaching the value y/l + ~ 1.27 and reaching 1 both at r — and r — > oo such that 



\/l + VT+4^ 
V2 



renders d(r) well defined at r — 1. Indeed, one can show that 



d{r) = 3 - 



V2 



111 1 



lnA(r) 



(21) 



(22) 



G 



In short, the prescribed A-shaped function A(r) absorbs the discrepancy of the function d(r) at r = 1 in order to make 
the d function well-defined everywhere including r = 1. For simplicity, one will choose 



A(r) - 1 + [^±2 - l]VFexp[l - yft. 



(23) 



One can easily show that this function peaks at r = 1 reaching the value \/l + \/E/\/2 and reaches 1 both at r = 
and r — > oo as required. And it has little effect elsewhere. This is exactly the function shown in Fig. 6. In addition, 
Figure 7 shows A(r) and g(r) together for comparison. 
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FIG. 6. A-shaped smooth function one needs is plotted as function r. Note that g and r are dimensionless quantities. 
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FIG. 7. A and g(r) are plotted as thick solid line, dashed line respectively. 



Therefore, one can plot the function of effective dimension d(r) in Fig. 8 
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FIG. 8. effective dimension d(r) is plotted as function r. 



Note that the effective dimension function d(r) does not represent the physical dimension of the theory. It only 
tells us roughly and effectively how the gravitational field g varies as r increases. In one aspect, one could view the 
potential as distortion of the geometry and the effective potential as the stiffness of the distorted geometry. Since 
the 1/r effect holds in all direction at large scale, not just along the disk plane, this behavior is quite different from 
simple geometric and dimensional deformation. All evidences indicate this sort of geometry deserves more attention. 



In this paper, one proposes a different view of MOND by looking at the physics related by the mass distribution 
and the coordinate scale of the system. In addition, one finds it useful to separate the effective gravitation field g(r) 
into a small scale (or short-distance )g s field and a large scale (or a long-distance) gi field that should be helpful for 
a better understanding of the underlying physics. The effective potential is obtained accordingly. 

We have also studied the changing pattern of the effective dimension d(r) which will be defined as a function of 
the physical scale r. One plots this d(r) functions hoping for a better understanding of the changing pattern of the 
effective dimension. Possible speculation with the Kaluza-Klein theory is also discussed in this paper too. 
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